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ABSTRACT 

We discuss new observational constraints on the abundance of faint high rcdshift Lyman-a emitters 
secured from a deep Keck near-infrared spectroscopic survey which utilizes the strong magnification 
provided by lensing galaxy clusters at intermediate redshift (z = 0.2 — 0.5). In each of 9 clusters we 
have undertaken a systematic 'blind' search for line emission with NIRSPEC in the J-band within 
carefully-selected regions which offer very high magnifications (> x 10-50) for background sources 
with redshifts z ~ 10. The high magnification enables the detection of emission at unprecedented 
flux limits (10 41 - 10 42 erg s _1 ), much fainter than those of conventional narrow band imaging and 
other spectroscopic searches. As the comoving volumes probed are small, our survey is designed to 
address the important question of whether intrinsically feeble star forming galaxies could provide the 
dominant ionizing flux at z ~ 10. Our survey has yielded six promising (> 5er) candidate Lyman-a 
emitters which lie between z — 8.7 and z = 10.2. We carefully discuss the validity of our detections 
in the light of earlier, apparently false, claims. We then focus on the likelihood that the detected 
line is Lyman-a from a distant lensed-source. All but one of the candidates remain undetected in 
deep HST optical images and lower rcdshift line interpretations can be excluded, with reasonable 
assumptions, through the non-detection of secondary emission in further spectroscopy undertaken 
with LRIS and NIRSPEC. We demonstrate that rigorous confirmation of the Lyman-a hypothesis 
for all of our candidates will be a challenge with current facilities. Nonetheless, as a result of our 
tests, we argue that at least two of our candidates are likely to be at z ~ 10. Given the small survey 
volume, this suggests there is a large abundance of low luminosity star-forming sources occupying 
the era when the reionization process is thought, from independent evidence, to be ongoing. While 
the predicted reionization photon budget depends upon a large number of physical assumptions, our 
first glimpse at the z ~ 10 Universe suggests that low luminosity star-forming galaxies contribute a 
significant proportion of the UV photons necessary for cosmic reionization. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
high-redshift — gravitational lensing 



1. INTRODUCTION 

After the epoch of recombination at z~1000, neutral 
hydrogen permeated the universe. The earliest observ- 
able structure dates back to the beginning of this era at 
the surface of last scattering traced by the cosmic mi- 
crowave background (CMB). The next most distant ob- 
servable represent the populations of ga laxies, gamma- 
ray bursts, and quasars at z ~ 6.2 — 7.0 (iHu et al. 2002 : 



Ive et al.1l200llTotani etafl 2006; K ashikawa et al.ll2006 



Fan et al.ll2006t iBouwens fc Illingworthl 120061 ). The in- 
terval between z ~ 1000 and z ~ 6.5 contains many 
landmark events: the formation of the first stars, the as- 
sembly of first galaxies, the growth of the first supermas- 
sive black holes that power quasars, and the reionization 
of neutral hydrogen in the intergalactic medium. One 
of the primary driving forces in observational cosmology 
today is to illuminate the growth of these first structures 
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in this important, but poorly understood, era. 

Recently, important constraints have emerged on the 
extent of star formation activity that occurred in the 
latter part of this era. The WMAP temperature- 
polariz ation cross correlati on signal on large angular 
scales (|Spergel et a l. 2006) implies scattering of mi- 
crowave photons by free electrons from ionizing sources 
at 2 ~ 10 — 12, suggesting an early period of effi- 
cient star formation. Meanwhile, Hubble Space Tele- 
scope (HST) and Spitzer observations have uncovered 
a population of very massive (M ste ;/ a r ~ 10 10 - 10 11 
M Q ) galaxies at z ~6 with well-establ i shed (>100 Myr) 



stellar populations (lEgami et al.ll2005t lEvles et al 



Mobasher et al. 2005; Yan et al. 2005; S tark et al 



2005 



2006 



Evles et all 120061: lYan et al.1 12006( 1. The assembly of 



such massive galaxies by these redshifts requires signif- 
icant star formation at yet earlier times. Observations 
of z-dropouts and J-dropouts with HST indicate a de- 
cline in the star formation activity of luminous Lyman- 
break galaxies out to z ~ 7 — 10 (jBouwens et al.l [20051 : 
IBouwens fc Illingworthl 2006): if there is significant star 
formation activity in this redshift interval, these observa- 
tions indicate it likely occurred in lower-luminosity sys- 
tems. 

Despite these indications of early activity, current fa- 
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cilities are poorly-equipped to undertake conventional 
searches for star-forming sources beyond z ~ 7. Even 
with HST and 8 meter class telescopes, continuum 
('drop-out') and narrow-band imaging Lyman-a searches 
probe only the most luminous, rare, sources at these 
redshifts; the contribution to the reionization flux from 
these sources may well be small. While narrowband sur- 
veys have had significant success a t dete c ting Lyman- 
q emitters at z ^ : 6 (iHu et all 119981 120021 12004 
Malhotra fe Rhoadsl [2004 IHu et alj 120051: I Stern et all 
20051: iKashikawa et alj|2006t IShimasaku et al.ll2006h only 
one Ly man-a emitter h as been convincingly detected at 
z ~ 7 (jive et al.l 12006). At z ~ 9, narrowband sur- 
veys for relatively luminous Lyman-a emitters (> 10 42 ' 5 

erg s" 1 ) have yet to fin d any sources ([Willis fc Courbinl 
120051 ICuby et al.ll2006h . While a few candidate Lyman- 
break galaxies beyo nd z ~ 7 have been s u ggested from 
photometric data (jBo uwcns et al. 2005; Labb e~et al.l 
20061: iRichard et alJ[2006t IBouwens fc Illingworthl 120061 : 
Henry et al.l 120071 ) . none has yet been spectroscopically 
confirmed, due largely to the difficulties inherent to ob- 
serving in the near-infrared. 

Recent observations of the faint en d slope of the star- 
forming luminosity functio n at z ~ 6 ([Yan fc Windhorst 
12004 IBouwens et' all 12006) support the suggestion that 
the bulk of the integrated star formation at high redshift 
may arise fro m very low lumino s ity sou r ces, as expected 
theor etically (|Barkana fc Loebl 120011 : IWvithe fc Loebl 
2006). However, forming stars in low mass halos is com- 
plicated by a number of feedback effects, any of which 
could significantly decrease the star formation efficiency. 
Clarifying whether low luminosity systems are abundant 
at early times requires observations probing well below 
the sensitivity limits obtainable with current facilities us- 
ing conventional methods. 

The study of lower luminosity sources rendered visible 
by the strong gravitational lensing induced by foreground 
clusters is the only current means to advance this quest 
before the next generation of telescopes is available. Al- 
though a typical lensing cluster magnifies sources by x 5- 
10 over 1 arcmin 2 for sources at z > 7, faint signals can be 
boosted > x 15 on the so-called critical line - the location 
of the lensing caustic in the image plane. Spectroscopy 
and detailed studies become feasible at otherwise impos- 
sible (unlensed) limits. F or several object s recently lo- 
cated via this technique (lEllis et al.l I2Q0TI: iKneib et al.l 
12004 ISantos et al.ll2004 lEgami et al.ll2005h . stellar con- 
tinuum slopes, Lyman-a profiles and star formation rates 
have been determined for sources whose intrinsic flux is 
close to th e faintest direct imagin g limits so far reached, 
I AB ~ 30 (|Beckwith et all [2006:. Until the availability 
of more powerful ground and space-based facilities, scan- 
ning the critical lines of massive clusters provides possi- 
bly the only reliable means of estimating the abundance 
of low luminosity star forming sources in the 7< z <12 
interval over which cosmic reionization is thought to take 
place. 

This paper presents the results of a spectroscopic sur- 
vey for strongly-lensed Lyman-a emitters at 8.5<z<10.4 
using the NIRSPEC near-infrared spectrograph on the 
Keck II telescope. The survey follows logically from a 
very similar, successful one undertaken at z ~ 4 — 6 us- 
ing the LRIS optical spectrograph on the Keck I telescope 



(|Santos et alj |2004) to which the reader is referred. 

A plan of the paper follows. We contrast our lensing 
search technique with other, more conventional probes of 
early star forming sources in §2. The NIRSPEC obser- 
vations and reduction methods are presented in §3. Here 
we describe the candidate selection criteria and derive the 
survey sensitivity function. A set of 6 promising candi- 
date z ~ 10 Lyman-a emitters arc presented in §4. Here 
we discuss in detail both the reality of the line detections 
and the likelihood that the line identifications represent 
high redshift Lyman-a. Noting the likelihood that a sig- 
nificant fraction of our sources may lie at z ~ 10, we in- 
fer the abundance of high redshift low luminosity sources 
and compare these with other constraints in §5. In §6 we 
discuss whether the abundance so derived might repre- 
sent a significant fraction of that necessary to complete 
cosmic reionization. We present our conclusions in §7. 

Throughout the paper, we have assumed a flat universe 
and (fi m , Qa)=(0.27,0.73) followin g the results pr esented 
in the initial WMAP data release ([Spergel et~afll2003h . 

2. GRAVITATIONAL LENSING AND SURVEY STRATEGIES 
- A CRITIQUE 

Strong gravitational lensing by intermediate redshift 
galaxy clusters offers the capability of extending the cur- 
rent generation of high redshift galaxy surveys so that 
lower luminosity sources can be brought into view. This 
added boost in sensitivity is particularly effective at the 
highest source redshifts (z > 7) where the bulk of star 
formation act ivity may be taking pla ce in low mass dark 
matter halos ([Barkana fc Loeb1l2001h . 

Lensing increases the apparent area of a background 
source while conserving the surface brightness; thus the 
gain in sensitivity is the magnification factor, M.. Rich 
clusters of galaxies at intermediate redshifts (e.g. z ~ 
0.2 — 0.5) amplify high redshift sources by factors of typ- 
ically M. ~ 3 — 5 over areas of a few arcmin 2 . However, 
boosts of M. ~ 10 — 50 are possible near the critical line 
- the theoretical locus of points in the image plane cor- 
responding to the caustic of infinite magnification. Clus- 
ters are usually characterized by both an inner and an 
outer critical line, the location of which varies as a func- 
tion of source (background) redshift. 

The accurate location on the sky of these critical lines 
can only be determined in the case of clusters for which 
detailed mass models are available. Such a mass model 
is primarily constr ained by the posit ions and redshifts 
of lensed features (jKneib et al.l I1996T ). A mass model 
can be validated, for example, by securing spectrosopy 
of lensed features whose redshifts were pre dicted by the 
mass model using an 'inversion' technique ([Ebbels et al.l 
1998). Currently, such high quality mass models are only 
available for ~20 clusters; however t his should increase i n 
the future with additional surveys (lEbeling et alJl200l . 

However, the gain in sensitivity afforded by strong lens- 
ing comes with a loss in survey area, compared to an 
equivalent unlensed survey, by a factor of M . This is an 
important point as the area of high magnification (say 
A4 >10) is already quite small (~ 0.1 arcmin 2 ). Pro- 
vided, as seems reasonable (§1), that the slope of the 
faint-end slope of the luminosity function is very steep, 
then strong lensing surveys should be a very effective way 
to characterize the amount of star formation at z > 7. 

Observers are considering several ways to conduct sur- 
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veys for z >7 Lyman-a emitters. Before describing the 
present survey, we therefore consider it instructive to ex- 
amine the benefits and drawbacks of the different strate- 
gies. 

The primary surve y techniques u s ed th us far are 
narrowba nd imaging fail et all l200l [2004 and spec- 
troscopy (|Martin et al. 120061 ). Neither technique has yet 
been routinely used in the near-infrared spectral region 
appropriate for z >7 and both methods can, in prin- 
ciple, be combined with strong lensing by focusing on 
fields containing rich clusters. At optical wavelengths, 
unlensed narrow-band imaging techniques typically reach 
flux limits of 2xl0~ 17 ergs cm -2 s , whereas spec- 
troscopic technique s probe to 5xl0 -18 ergs cm -2 s _1 
(jSantos et al.ll2004h . 

In the near-IR, a longslit spectrograph such as NIR- 
SPEC offers a ^8 A resolution over a moderate redshift 
interval (Az=2). Although bright atmospheric emission 
lines restrict visibility for 50% of the redshift range in the 
J-band, for a resolving power R~1500, the gain in red- 
shift coverage is still a factor of 10-100 over the narrow 
band technique (Az—0. 01-0.1). Per exposure, there is 
however a very small field of view (~30 arcsec 2 ) c.f. nar- 
rowband imaging (~ 10 arcmin 2 for a cluster but much 
larger for panoramic field surveys). 

Assuming the characteristics for current instrumenta- 
tion and their detectors, we compare the efficiencies of 
the two methods in the lensing case in Figure 1. Details 
of the computation are given in the Appendix. Both 
lensed survey techniques probe at least an order of mag- 
nitude deeper than conventional near-IR Lyman-a nar- 
rowband surveys. Of the two lensing methods, long-slit 
spectroscopy is better suited for detecting the faintest 
Lyman-a emitters for detectors 1024 x 1024 pixels in size. 
This can be explained as follows. While a narrowband 
imager covers a larger area on the sky, only a very small 
fraction of this area is highly magnified (e.g. >10x). The 
small gain in the highly-magnified survey area is insuffi- 
cient to compensate for the significantly smaller redshift 
coverage. Moreover, the ~100 A bandwidth typical for 
many narrowband filters results in significantly poorer 
sensitivity limits. Taking the various factors into consid- 
eration, a long slit spectroscopic survey represents the 
more efficient technique for detecting sources with star 
formation rates below 10 42 erg s _1 , and hence is the op- 
timal strategy for determining whether reionization was 
predominantly caused by large numbers of feeble sources. 

3. NIRSPEC CRITICAL LINE SURVEY 

3.1. Cluster Sample 

Our goal is to constrain the abundance of sources up 
to an order of magnitude below current survey sensitiv- 
ity limits. Accordingly, following the discussion above, 
we have conducted a spectroscopic survey for Lyman- 
a emission along the critical lines of 9 clusters at z ~ 
0.2 — 0.5. We limited our observations to those clusters 
with mass models that are well-constrained by the com- 
bination of HST imaging and ground-based spectroscopy 
of multiply imaged systems. The list of survey clusters 
and the associated mass models are summarized in Table 
1. 

Many criteria entered into the final selection of our 
9 clusters. An essential criteria was the availability of 
deep HST imaging for constructing samples of multiply- 



imaged systems (|Kneib et al.l Il996t iBroadhurst et al.l 
2005). The most well-understood clusters in our sam- 
ple contain 10-30 multiply- imaged systems, a signi- 
fication fraction of which have spectroscopic redshift 
(e.g. Abell 1689 has 33 such systems of which 21 are 
spectroscopically-confirmed, Richard et al, in prep). This 
data is necessary to astrometrically pinpoint the location 
of the critical line along which the NIRSPEC slit is ori- 
ented. 

Deep HST optical imaging is also essential for verify- 
ing, via the absence of any detectable continuum, that 
any lensed emitters are genuinely at high redshift. Due 
to neutral hydrogen absorption we would not expect any 
of our z ~ 10 emitters to be seen at optical wavelengths 
and HST provides the deepest verification. 

A final consideration is the availability of deep opti- 
cal spectroscopy which is helpful in order to check for 
associated emission lines that might arise if any NIR- 
SPEC line detection arises from a lower redshift source. 
Giv en the comprehensi ve critical line survey undertaken 
by (jSantos et a l. 2004) using LRIS, where possible we 
chose the same clusters and explored the same critical 
line regions. 

In practice, telescope scheduling and weather losses 
compromised a perfect adoption of the latter criteria. 
However, for each of the 9 clusters, the mass model, the 
critical line location and associated magnification prop- 
erties are well understood. 

3.2. Observations and Data Reduction 

We utilized the Near I nfrared Spectrometer (NIR- 
SPEC, [McLimirS] [TUl) mounted on the Keck II 10 
m telescope at Mauna Kea to perform our survey. Obser- 
vations were performed in the J-band (1.143-1.375 /im) 
with a slit 42" long and 0'.'76 wide. We used the low- 
resolution mode of NIRSPEC with the 75 line mm -1 
grating which offers a net resolving power of R~1500 
and a spectral resolution of ~8 A. At the wavelengths 
sampled within the NIRSPEC J-band filter, any detected 
Lyman-a emission would correspond to sources with red- 
shifts between z = 8.5 and z = 10.4. 

Figure 2 shows how the targeted slit positions match 
the predicted location of the z ~9 critical lines for each 
of the 9 clusters in our survey. To maximize efficiency, 
where a choice was available we selected regions where 
the critical line is straight enough to provide high mag- 
nification across the entire slit. 

Typically, 2-5 slit positions were observed per cluster. 
For each slit position, we took 6-10 ten minute exposures 
using a three point dither pattern in which the telescope 
was offset ±3" along the slit. Table 2 summarizes the 
details of our NIRSPEC observations. In four observing 
runs, we observed 35 slit positions corresponding to a 
total sky area of 0.3 arcmin 2 . 

Spectra were flat-fielded and sky-subtracted using IDL 
scripts written by G. Becker (2004, priv ate communica - 
tion). Following techniques described in Kelson (2003), 
the two-dimensional spectra were not straightened prior 
to sky subtraction, thus ensuring sub-pixel sampling of 
atmospheric emission lines which, in turn, led to a sig- 
nificantly improved removal of the sky background. 

The camera distortion and spectral curvature were 
computed by fitting standard star traces along the slit 
and atmospheric emission lines across the dispersion axis. 
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Knowledge of the distortion and spectral curvature was 
used to construct arrays that provide proxies to the slit 
position and wavelength for each exposed pixel on the 
NIRSPEC detector. Using this information, sky back- 
ground modeling was performed using a two-dimensional 
b-spline fit, where a low-order polynomial is fit to the slit 
illumination and a b-spline is fit along the dispersion axis. 
The wavelength calibration is subsequently computed us- 
ing atmospheric emission lines. For each exposure, we 
also compute 2D-variance arrays by summing the contri- 
bution from the dark current, flat-field, read noise, sky 
background, and source counts to the variance. 

Following techniques introduced by Santos et al (2004), 
we determined the astrometric position of each spectral 
exposure by registering images from the NIRSPEC slit- 
viewing camera (SCAM) to HST images of the same 
field. WCS coordinates were computed for each exposed 
pixel on the detector using the slit position grid described 
above. Since these are galaxy cluster fields, many bright 
objects appear in the 46x46 arcsec 2 field of view, en- 
abling accurate registration on to the WCS of the HST 
images with a typical rms of 0'.'15 (less than 25% of the 
slit width). Offsets between exposures of a given slit po- 
sition are calculated from the registered SCAM images, 
and the 2D-spectra and their associated 2D-variance ar- 
rays are subsequently shifted and combined. To remove 
cosmic rays and bad pixels, we median combine the data 
rejecting the brightest and faintest frame at each pixel. 

Our ability to detect faint emission lines is strongly 
dependent upon the accuracy with which the offsets be- 
tween the different exposures of a slit position are known. 
If a bright emission line from a foreground galaxy lies 
serendipitously on the slit, the accuracy of the offsets can 
be determined by comparing the shifts derived from the 
SCAM images to the offsets between the centroid of the 
emission line. Suitable emission lines are present in sev- 
eral of our slit positions and applying the aforementioned 
test, we find the offsets determined from the SCAM im- 
ages are typically good to Of.'l. At less than 15% of the 
slit width, this uncertainty does not compromise the S/N 
of faint emission features. 

One additional concern in dealing with faint emission 
line objects is the possible transverse drift of the object 
across the slit from frame to frame when the seeing is 
less t han a slit width (|Pell6 et alj|2004 : IWeatherlev et all 
20041). This would cause the centroid of the emission line 
to drift in the dispersion direction of the spectrun. Un- 
corrected, this would increase the area over which the 
emission line flux is spread, thereby reducing the ob- 
served S /N of the emission line in the final stacked spec- 
trum. We quantify the magnitude of the transverse drift 
by measuring the frame-to- frame centroid of a very bright 
[OIII] emission line in NIRSPEC observations taken with 
exceptional seeing (0'.'45). We find that the standard de- 
viation in the centroid position over 9 frames is 0.8 pixels, 
or 0'.'2. Adding this measured tranverse drift in quadra- 
ture to the seeing FWHM only increases the FWHM of 
the line by 5%. Given that the seeing was never better 
than the value quoted above, it is clear that the drift has 
a negligible effect on the detectability of emission lines. 

3.3. Survey Sensitivity 

We now determine our survey sensitivity function. We 
first discuss the limiting line flux as a function of wave- 



length and source redshift. Converting this flux to a lim- 
iting source luminosity requires knowledge of the mag- 
nification across the NIRSPEC slit. In practice, we will 
quote our limit in terms of that line flux we can expect 
to detect, at 5ct confidence, for an unresolved emitter 
whose rest-frame line width i s typical of a fain t, but very 
well-studied, z ~5.7 system (jEllis et alj|2001l) . 

To accomplish this, we computed the total noise as a 
function of dispersion and slit position using the vari- 
ance arrays obtained in the data reduction (§3.1). This 
comprises the sum of the variance from the sky back- 
ground, read noise, flat-fielding, dark count as well as 
the source counts. Flux calibration of the variance ar- 
rays was performed using observations of spectroscopic 
standard stars. The 5cr limiting Lyman-a line flux was 
then computed as a function of slit position and wave- 
length by calculating the signal that is 5 times the root 
of the variance in an aperture whose spatial dimension is 
twice the measured seeing disk (typically 0'.'7) and whose 
spectral dimension equals the Lyman-a line width, ~300 
k m s -1 of the lensed system studied with high precision 
bv lEllis et al] (l200lh . 

We find the median 5<r limiting line flux across a typi- 
cal slit position with 1.5 hours of integration is 2x 10~ 17 
erg cm -2 s _1 (see Figure 3); the sensitivity varies be- 
tween 0.9-3 x 10~ 17 erg cm" 2 s" 1 depending on the 
proximity to atmospheric OH lines. This flux limit is 
similar to that reached in the z = 8.8 narrow band survey 
of GOODS-S discussed in rCubv et alJ (12001 (1.3xl0~ 17 
erg cm -2 s _1 ) and slightly less sensitive than the z = 8.8 
narrowband survey presented in lWillis fc Courbml (|2005f ) 
(3.3xl0~ 18 erg cm -2 s _1 ); however the lensing magni- 
fication enables the detection of sources that are signifi- 
cantly less luminous than those detectable in either con- 
ventional narrowband survey. 

3.4. Cluster Magnification 

As discussed i n our precursor optical paper 
([Santos et al.l l2004h . the cluster mass models listed 
in Table 1 provide a 2-D map of magnifications for 
any given source redshift. The details by which these 
magnification s are realized foll ow techniques described 
in detail by (jKneib et al.l fl996). Further discussion of 
the code, LENSTOOL, used to create the mass models 
for systems for which rich d e tail i s available Julio et 
al, (in prep and ISmith et al.l ((2005). From each mass 
model, the magnification can be computed as a function 
of source redshift and position. 

In practice, for each position and wavelength, the mass 
models produce a matrix of the convergence and shear. 
The magnification is then determined from the conver- 
gence k and shear 7 using the following equation: 

^ (n ' z)= [i-K(n,z)P-7(n,^ (1) 

where CI is the position on the sky and z is the source 
redshift. 

A key question that often arises in the consideration 
of critical line surveys is the accuracy with which the 
location of the line can be pinpointed and also that of 
the associated magnification map. To address this ques- 
tion, we have used a Markov chain Monte Carlo sampling 
algorithm within LENSTOOL (Julio et al., in prep) to 
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map the probability density in the mass model parameter 
space for each cluster. We then generated a magnifica- 
tion map for each model and in turn computed the mode 
of the distribution at each pixel as well as the range of 
magnification factors which bracket the mode by ±34%. 
This we take as the uncertainty in the magnification. The 
errors derived for each cluster are applied in the deter- 
mination of both the candidate luminosities (§5) and the 
survey volume (§6.1). 

Taking each exposed pixel on the NIRSPEC detector 
as an element of solid angle, the convergence and shear is 
computed as a function of position via interpolation be- 
tween the NIRSPEC coordinate grid and that of the con- 
vergence and shear maps. For each NIRSPEC pixel, the 
associated redshift at which a Lyman-a emitter would 
be located is calculated using the wavelength solution 
described above. The convergence, shear and Lyman-a 
redshift is then inserted into equation 1 and the mag- 
nification M. at each pixel determined. Likewise, the 
magnification uncertainty map is interpolated onto the 
NIRSPEC coordinate grid enabling the computation of 
the error as a function of slit position. Both the magni- 
fication and associated error maps were smoothed using 
a l'.'O median box filter to ensure that our results are 
not corrupted by sharp high peaks in the magnification 
distribution. The change in the magnification error with 
redshift is negligible; hence we adopt the uncertainity 
at z = 9 to be the magnification error for all Lyman-a 
redshifts probed by our survey. 

As an example of the magnifications sampled, the 
wavelength-independent magnification along one of our 
survey slit positions is displayed in Figure 4 for three dif- 
ferent redshifts. The magnification is > x20 over nearly 
the entire slit and fairly constant with redshift. We also 
plot the associated la uncertainty as a function of posi- 
tion along the slit at one of the redshifts. Even allowing 
for the uncertainty in the lensing model, the magnifi- 
cation is still expected to be uniformly high across the 
entire slit. This is largely a result of the fact that we only 
select clusters with well-defined mass models. Further- 
more, in regions where the location of the critical line is 
known less precisely, there is generally a region ~2'.' in 
width around the critical line where the magnification 
factor is greater than fifteen, so the positioning of the 
NIRSPEC slit is always in a region of high magnifica- 
tion. 

Figure 5 presents a histogram of the magnification in 
our survey. Throughout the entire survey area, the me- 
dian magnification provided to background sources is 
~ x20, a factor o f two larger t han in the optical sur- 
vey conducted by ISantos et al.l (|2004f ) . The increased 
magnification is due to several factors. First, ultra-high 
magnification regions (x 15-30 throughout the slit-area) 
were preferably selected to allow very low-luminosity sys- 
tems (< lO^ergs^ 1 ) to be probed. Second, since the 
NIRSPEC longslit covers a factor of three smaller an- 
gular extent on th e sky than the LRIS longslit used in 
ISantos e~a l. (2004), we were able to avoid regions where 
the critical line curved off the slit, resulting in high mag- 
nification uniformly across the slit. 

3.5. Limiting Lyman-a Luminosity 

Given our knowledge of the limiting Lyman-a flux 
(§3.2) and the spatially-dependent magnification (§3.3), 



we are now in a position to calculate the limiting source 
Lyman-a luminosity, Ln m (f2, z). This we define as the 
least luminous source detectable (with 5c confidence) in 
a given volume element, allowing for the cluster magni- 
fication M.. 

Allowing each NIRSPEC pixel to correspond to a vol- 
ume element, we can calculate the limiting luminosity 
across each slit in the survey according to the following: 



£n m (ri, z) = 



^(\ + zfDl{z) fa m (z) 



T(n) 



M(n,z) 



(2) 



where D c (z) is the comoving distance of the volume ele- 
ment and T(J7) is the slit transmission. 

In the absence of the lensing magnification, the typical 
limiting Lyman-a luminosity is 2x 10 erg s _1 in be- 
tween atmospheric OH lines. However, more than 50% 
of the survey area is magnified by > x30, enabling the 
detectability of sources as faint as 7x 10 41 erg s" 1 over 
~0.15 arcmin 2 in the image plane. Conventional narrow- 
band Lyman-a surveys at z > 6 typically reach limiting 
luminosities between 10 42 and 10 43 erg s _1 . Our NIR- 
SPEC survey probes over an order of magnitude fainter 
and thus can clearly provide constraints on star form- 
ing sources that are otherwise out of reach of traditional 
methods, at least with current facilities. 

4. IDENTIFICATION AND REALITY OF THE CANDIDATE 
LYMAN-a EMITTERS 

At the faint limits now being probed, we have found the 
reliable identification and verification of distant Lyman- 
a emitters to be a very challenging endeavor, even with 
the most powerful facilities available to us. The bur- 
den of proof that a detected line is (a) real and (b) 
truly arises from a highly-rcdshifted Lyman -a emission 
line is very great given earlier controversies (|Pell6 et al.l 
l200i iBremer et al.ll200i iSmith et al.1[200l . So far as 
the identification process is concerned, the key issue is 
to ensure candidate emission features are not spurious 
and do not arise from detector artefacts lWeatherlev et ail 
(2004). Concerning the redshift verification, a new set of 
challenges emerge at z ~10 since additional (confirma- 
tory) emission lines are not available from ground-based 
facilities (see §5). 

Beginning with the identification process, each of the 
sky-subtracted 2-D spectra were independently inspected 
for Lyman-a candidates by three of the authors (D.P.S, 
R.S.E, and J.R.). We demanded that any candidate fea- 
ture be both extended over multiple pixels (in order to 
differentiate a signal from bad pixels and other detec- 
tor artifacts) and located away from contaminating OH 
airglow lines. The three lists of line candidates each con- 
tained ~ 25 possible features with many being in com- 
mon between two of the three lists. A final catalog of six 
promising candidates was constructed by selecting those 
features with the highest significance common to all three 
observers. These are displayed in Figure 6. The S/N ra- 
tios of these candidate emission lines, defined as in §3.2, 
varies between 5 and 8 (see Table 3). 

A spurious feature could arise if the reduction process 
does not adequately remove a cosmic ray or bad pixel. 
Obviously dithering the telescope so the candidates move 
up and down the slit eliminates most detector anomalies. 
However artefacts can remain if the spectra are rebinncd 



G 



Stark et al. 



during the wavelength calibration; in this case, such bad 
pixels woul d be smoothed out and a ppear similar to emis- 
sion lines ([Weather lev et al.ll2004fl . Fortunately, we do 
not face this problem in our Keck dataset since our re- 
duction technique altogether avoids rebinning the spec- 
tra (see §3). Nevertheless, one may still worry that the 
measured emission line flux originates from one or two 
frames with overlapping cosmic rays or bad pixels. We 
have tested this possibility by median combination of all 
of the exposures of slit positions containing our candi- 
dates and rejecting the three highest pixels. We also me- 
dian combined half of the exposures (chosen randomly) 
for each candidate. In both tests, all six candidates are 
still clearly visible (albeit at lower S/N for the latter 
test), suggesting that our emission features do not arise 
from flux in only one or two frames. 

Finally, it is conceivable that the noise distribution of 
the spectral data is such that 5cr noise features are much 
more common than would be expected under Gaussian 
statistics. Such a curious situation might arise if there 
was some unforeseen property of the NIRSPEC detector. 
In this case our candidate emission features might simply 
be peaks in the noise distribution. We investigated this 
possibility by computing the flux in apertures centered at 
random positions in areas of the spectra where the sky- 
subtraction is clean. The aperture size was matched to 
the expected size of a. Lyman-a feature at high redshift: 
the spatial FWHM was taken to be the size of a typical 
seeing disk, 0'.'6 and spectral FWHM corresponding to 
typical line width of Lyman-a emitters, 300 km s _1 ) 

A normalized histogram of the distribution of summed 
aperture fluxes is presented in Figure 7. The frequency 
with which a feature with S /N = 5 appears in our spectra 
is ~0.02%, which corresponds to just under a 4cr detec- 
tion for a Gaussian distribution. Though the frequency 
of features with fluxes comparable to our candidates may 
be slightly higher than indicated by Gaussian statistics, 
the test clearly illustrates that features of prominence 
similar to those of our candidates are extremely rare in 
clean regions of the spectra. 

5. VERIFYING THE CANDIDATE REDSHIFTS 

If the emission features discussed above are Lyman-a 
, then the source redshifts for our 6 candidates range be- 
tween z — 8.7 and z = 10.2 and the observed fluxes 
span 2-5xl0~ 17 erg cm" 2 s -1 . Converting these ob- 
served fluxes into intrinsic luminosities is complicated by 
uncertainties in the fraction of emitted Lyman-a photons 
that reach the observer. We can compute a lower limit 
to the unlensed Lyman-a luminosity by converting the 
fluxes to luminosities and dividing by the magnification 
factor. With these assumptions, the derived unlensed 
luminosities range from 1.6-50xl0 41 erg s -1 . 

We convert these Lyman-a luminosities to star for- 
mation rates assuming case B recombination: L Q = 
|hi/Lya(l — fcsc)Q(H)M^ where hi/ a is the energy of a 
Lyman-a photon, Q(H) is the Hydrogen ionizing pho- 
tons flux, and f esc is the fraction of ionizing photons 
that escape from the galaxy (and hence are not emitted 
in recombination emission lines). Assuming a Salpeter 
initial mass function, l/20th solar metallicity, and a min- 
imum and maximum stellar mass of 1 and 20 M , respec- 
tively, the hydr ogen ionizing photon flux is 4xl0 53 s _1 
(|Schaererll2003f) for galaxy forming stars continuously at 



a rat e of 1 M Q yr _1 . For f esc <C 1 (e.g. IShaplev et al.1 
2006), this results in a conversion factor of 1 M Q yr _1 
= 4.2xl0 42 erg s" 1 . With these assumptions, the star 
formation rates range from 0.038 to 1.2 M Q yr _1 . If a 
Scalo IMF was adopted instead, th e star formation r ates 
would be a factor of three greater (|Loeb et al.ll2005[ ). 

We note that these star formation rates are uncertain 
for several reasons. If dust in the galaxy absorbs any 
Lyman-a photons, then the intrinsic luminosities and 
corresponding star formation rates will be larger. Like- 
wise, if only a fraction of the Lyman-a line is transmitted 
through the IGM due to the presence of neutral hydro- 
gen, then Lyman-a luminosity and star formation rates 
will again be underestimated. On the other hand, if the 
stellar IMF at these early times is more top-heavy than 
a Salpeter or Scalo IMF, more ionizing photons are pre- 
dicted per unit of star formation. Similarly, if the metal- 
licity of these galaxies is less than what is assumed above, 
the efficiency of ionizing photon pro duction (and hence 
Lyman-q luminosity) is increased ([Tumlinson fc ShulU 
I2000t iTumlinson et all 120011 : ISchaereij|2003[ ). Hence, in 
these cases, the star formation rates tabulated will be 
overestimates. 

But how can we be sure, with a single emission line, 
that we have truly detecteda Lyman-a emitter? Nor- 
mally , for a narrow-band imaging survey (Hu et al.l l2002t 
1 2004 iMalhotra fc Rhoadsl 12004 iShimasaku et al.1 120061 : 
iKashikawa et al.ll2006h or an unlensed spectroscopic sur- 
vev (|Martin et al.ll2006"f ) , contamination from foreground 
[O II] and Ha emitters can be significant. However, it 
is worth pointing out that such foreground contamina- 
tion is likely to be much reduced for a critical-line survey 
designed to probe faint emitters since, for lensing clus- 
ters at z ~0.2-0.5, all contaminating emitters will also 
be lensed and thus their respective most-likely location in 
the image plane will be spatially offset with respect to that 
for a source at z ~2. The degree to which this benefits a 
lensed survey with respect to a blank field survey will de- 
pend on the relative magnifications at each redshift and 
the shape of the faint end of the respective line emitting 
luminosity functions. 

A redshift that the general community normally ac- 
cepts is usually one in which at least two, preferably 
more, features are reliably identified. Single line identifi- 
cations are quite naturally questioned as unreliable. As 
we probe to z ~ 10 with ground-based telescopes, how- 
ever, we reach the situation where multiple-line identifi- 
cations are simply not practical. With Lyman-a in the 
J-band no additional feature is visible in the full practical 
spectroscopic range from the UV limit to the beginning 
of the thermally-dominated background at 2y«m. The In- 
frared Spectrograph (IRS) on Spitzer in principle offers 
the possiblity of detecting Ha emission from z ~ 10 ob- 
jects. However, sensitivity calculations (Lacy et al., in 
prep.) suggest that the Ha/Lyman-a flux ratio would 
have to be significantly greater than that expected for 
case B recombination (perhaps due to Lyman-a suppres- 
sion from resonant scattering and absorption) in order for 
this to be an effective verification. 

A second-best verification of Lyman-a that might be 
convincing would be the identification of an associated 
Lyman brea k in broadband ima ging of the continuum 
distribution (|Stanwav et alj [20041) or the presence of a n 
asymmetric line profile ([Ellis et aLlfeOOltlHu et"aiT 2002). 
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The Lyman break is due to intergalactic absorption of 
flux shortward of the Lyman-a line and causes all z > 7 
sources to be undetected at optical wavelengths. How- 
ever, the lack of an optical broad-band detection need 
not necessarily rule out a lower redshift source, for ex- 
ample those with very strong emission line spectra and 
weak continuum emission. 

Asymmetric line profiles are expected because of res- 
onant scattering by neutral hydrogen within the host 
galaxy and in the IGM. This leads to absorption on the 
blue side of the line as the light redshifts into resonance, 
leaving a sharp cutoff on the blue side of the line and 
an extended red damping wing. Unfortunately, the de- 
tection of an asymmetric line profile requires very long 
integrations and high spectral resolution in order to de- 
tect the damping wing. 

A key factor in our consideration of how to verify the 
Lyman-a label for our 6 candidates is the importance of 
exquisite conditions to see such faint emitters. The re- 
quired seeing and transparency occur only 40-50% of the 
time even on an excellent site such as Mauna Kea. Con- 
servatively, we estimated that to measure the Lyman-a 
profile of any one of our 6 candidates with adequate sig- 
nal/noise would take 8-10 hours of integration on Keck 
II. Allowing for weather, to follow-up all 6 candidates in 
this fashion would take 24 nights; clearly this is not a 
practical proposition! 

One might conclude we have reached new territory 
where we may never quite know with certainty whether 
an object is at z ~ 10. In practice, however, there is a 
way forward although it is less definitive and more statis- 
tical in nature. The approach we will adopt is to search 
spectroscopically in various passbands for lines associ- 
ated with alternative, lower redshift explanations for the 
detected J-band emission line. In ruling out a particular 
alternative identification, we must make some assump- 
tions about the likely emission line spectrum of the in- 
terloper. As such this is then a statistical process since 
we cannot cover all possible emission line ratios. While 
this method will never confirm the Lyman-a interpreta- 
tion with 100% certainty, it is clearly an essential pre- 
requisite to spending 2 nights per candidate measuring 
the line profiles. 

In the following subsections, we proceed to apply the 
tests discussed above to our candidates in order to con- 
strain their redshifts. In §5.1, we examine the stacked 
line profiles of the six candidates to see if we can collec- 
tively assess the likelihood the bulk represent Lyman-a 
emission; in §5.2, we test for the presence of a Lyman 
break in deep broadband photometry; finally in §5.3, we 
discuss our follow-up spectroscopic program to test for 
low-redshift interlopers. 

5.1. Stacked Line Profile of Candidates 

Asymmetric line profiles cannot be discerned in the in- 
dividual emission features, because of the low spectral 
signal/noise. As discussed, much longer integrations (8- 
10 hrs) would be needed to characterize the emission line 
profiles and this seems pointless until lower redshift in- 
terlopers have been rigorously tested. 

Nonetheless, one can ask whether statistically as a 
combined set, the combined line profile of all 6 candidates 
reveals an asymmetric profile suggestive of a Lyman-a 
origin. The practicality of such an exercise is limited 



by a number of factors. First the resolution of NIR- 
SPEC is only moderate; R -1500 c.f. R -15 000 f or the 
ESI spectrum obtained at z=5.7 (jEllis et al.ll2001h . Sec- 
ond, stacking assumes that all 6 lines are Lyman-a with 
similar profiles and, most importantly, requires an accu- 
rate registration using the central wavelength within each 
low signal/noise profile. Finally, as the adjacent OH sky 
lines are distributed differently for each emitter, wave- 
length regions that are clean in one candidate become 
averaged with noisy regions in another so the contrast 
of the stacked line does not increase as \/~N. Even with 
weighted addition, we did not find the test to be very 
conclusive. 

5.2. Searching for Lensed Pairs 

In several of the lensed emitters and drop-outs we 
have published in the past (jEllis et al.ll200lHKneib etahl 
l2004h . the location of a second or third counter-image 
has been a particularly convincing demonstration of both 
the lensed hypothesis and, via the mass model, the ap- 
proximate redshift. Indeed, the accurate recovery of the 
third image for the z ~ 6.8 drop-out source in Abell 2218 
(|Kneib et al.ll200l was one of the major arguments jus- 
tifying its location beyond z ~6. 

Generally speaking, image pairs are expected in strong 
lensing configurations depending on the degree of align- 
ment between the background source and the cluster 
caustics, the physical size of the source and, in this case, 
the ground-based seeing. The key feature which led to 
the satisfactory recovery of image pairs in our earlier 
studies was a broad-band detection with HST. As an 
example, in the case of the z = 5.7 pair discussed by 
Ellis et al (2001), only a single emitter was recovered by 
LRIS. The location of the second image required a detec- 
tion by HST. Indeed, pairs were not detected for any of 
the other lensed so urces in the subsequent LRIS survey 
(jSantos et al.ll2004h . 

Some true pairs may simply not be resolved in our 
ground-based spectra and, quite possibly, the second im- 
age lies outside our NIRSPEC slit coverage. For the most 
well-constrained mass model, the location of a counter- 
image can usually be predicted with an uncertainty of 
O'.'5-l'.', i.e. to within 1-2 NIRSPEC slit widths. Al- 
though it is practical to consider chasing the second im- 
ages, this could still be fairly time-consuming and we 
concluded that priority should be given in eliminating 
low-redshift interlopers from our catalog of candidates. 

5.3. Broadband Photometry of Candidates 

The next test is to determine whether optical contin- 
uum emission is visible for any of our candidates. Such 
a detection would clearly rule out a high z interpreta- 
tion although non-detec tion does not imply the opposite. 
Deep ACS/WFPC-2 (|Kneib et all l200l iSmith et all 
2002) images are available for all 9 clusters so we have 
examined these images at the precise locations of each of 
our 6 emitters. 

No definitive broadband optical detection is seen at the 
location of any of the candidates to R~27.5 (3<r, Fig- 
ure 8). One of the candidates (Abell 2219 c2) is very 
marginally detected in the z' FS5QLP -h&Tid; if this flux is 
associated with the J-band emission feature, then the 
high redshift interpretation would seem unlikely. Al- 
though 3 of the candidates (Abell 2219 cl, Abell 1689 
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cl, c2) are located within 2'.'6f brighter objects, in all 
cases, our registration sufficiently accurate to rule out 
the possibility that the J-band emission emanates from 
the bright objects. 

A detection in deep near-infrared imaging would 
be more interesting. Coupled with the optical non- 
detections discussed above, this might hint at the pres- 
ence of a Lyman break, lending credence to a z ~10 inter- 
pretation. In fact, none of the candidates is detected at 
or above the 5tr level (~27 AB mag) in deep broadband 
HST Jnow, Hi6ow and Keck J-band imaging. Abell 
2219 c2 is marginally detected at ~ 2cr in the Hi6ow- 
band. 

In fact the lack of definitive near-infrared detections 
for our Lyman- a candidates is not at all surprising 
if they are at high-z. If the equivalent width distri- 
bution of Lyman-a emitters at z ~ 10 is compara- 
ble to that at z ~4-6, then we would need to reach 
-ffi60W = 27.7 (1 magnitude below our sensitivity limit) 
to detect o ur candidates, if they are indeed Lyman- 
q em itters ()Malhotra fc Rhoads 2004; Shimasak u et al.l 
120011 . 

Thus, we can conclude that our deep imaging adds a 
further important component in the case for the hypoth- 
esis that bulk of our emission line sources lie at z ~10. 

5.4. Spectroscopically Testing Low Redshift Scenarios 

Aside from Lyman-a, the most likely alternative iden- 
tifications for the J-band emission features include Ha, 
[OIII], H/3, and [Oil], in which case the features would 
correspond to lensed sources at z = 0.5 — 2.5. Sources lo- 
cated at these redshifts would almost certainly have addi- 
tional, associated emission lines bright enough to be seen 
in the optical and near-infrared spectral regions (see Ta- 
ble 4). We have constrained each of these lower redshift 
hypotheses by undertaking further spectroscopy, placing 
stringent limits on the presence of lines. 

This technique has value so long as the expected line 
spectrum of the foreground interlopers is known. In prac- 
tice selective extinction and excitation differences make 
the relevant line ratios somewhat uncertain. The best we 
can do is to approach the problem statistically assuming 
a typical HII region spectrum (Figure 9) allowing some 
leeway in the line ratios. 

We generate the HII region spectrum by computing 
flux ratios of bright recombination lines (Lyman-a, Ha, 
and H/3) assuming Case B recombination and accounting 
for the possibility that a signficant fraction of Lyman-a 
photons are absorbed by dust thereby reducing its flux 
by up to a factor of three. With these assumptions, we 
adopt the following emission line ratios for hydrogen re- 
combination lines: Lyman-a : Ha : H/3 — 7.4-22 : 2.8 
: 1.0. Flux ratios involving the forbidden oxygen lines 
([Oil] and [OIII]) can be determined from observations of 
low- and high-metallicity galaxies (Pilyugin et al. 2000, 
2001). Here we adopt a ratio of [OIILU5007 : [OIII]A4959 
: [Oil] : H/3 = 0.8-7 : 0.3-2 : 0.8-3 : 1. Using these ratios 
and the conversion between Lyman-a luminosity and star 
formation rate (§5), we generate a template spectrum 
(Figure 9); the continuum flux spec trum is determined 
for a g iven star formation rate via a lBruzual fc Charlotl 
(2003) model with identical stellar properties to those 
assumed in §5. 

There are certainly a number of shortcomings in such 



an approach. Predicting the flux of Lyman-a is always 
difficult due to the resonant nature of the transition. 
While wc allow for a factor of three suppression in the 
Lyman-a flux, it is possible that a larger fraction of 
Lyman-a photons could be absorbed. Furthermore, the 
flux ratios, especially those involving the forbidden oxy- 
gen lines, are dependant on the metallicity and the effec- 
tive temperature of the ionizing stars. However, as we 
will show below, our observational flux limits are often 
sufficiently tight to rule out low-redshift hypotheses for 
a wide variety of metal abundance and stellar effective 
temperature assumptions. 

Following the above logic, after completion of the blind 
scanning survey, we set out to secure additional optical 
and H-band spectroscopy for all 6 candidates to apply 
this interloper-rejection method. Sadly, after two sea- 
sons of observing, weather and other observing vagaries 
mean that exhaustive consideration of all interloper pos- 
sibilities is only available for 3 of our candidates. For the 
other 3 candidates, only partial coverage is in place. 

Table 4 summarizes the various hypotheses for the J- 
band emission and lists, for each case, which additional 
lines would be expected and at what wavelengths. Op- 
tical spectroscopy is very efficient at constraining the 
likelihood that our J-band emission arises from either 
Ha (z~0.9) or [O III] (z ~1.5) since, in these cases, 
the LPJS data probes very deep in the associated wave- 
length regions where [O II] would be seen (at 7100 A 
and 9300 A respectively). In some cases, the survey J- 
band NIRSPEC data can also be used to constrain the 
hypothesis that the primary detection is H(3 or [O III] 
where we would expect to locate one or both of the [OIII] 
A4959,5007 pair and H/3. In practice, this is complicated 
by the fact that lines can be obscured by the OH night 
sky or seen at very low significance (e.g. 2a). We found 
that such cases can be more effectively dispensed with us- 
ing H-band spectroscopy to search for Ha which should 
be easily observable if it lies in between the sky lines. 

5.5. Best Candidates: Abell 68 cl and Abell 2219 cl 

The most convincing elimination of foreground inter- 
loper status has been achieved for the two candidates 
Abell 68 cl and Abell 2219 cl. In our subsequent analy- 
sis we will assume these emitters are Lyman-a and so we 
describe these spectroscopic constraints in some detail. 

In addition to the null broad-band detections discussed 
earlier, deep optical LRIS and H-band NIRSPEC spec- 
troscopy was secured in good conditions for both candi- 
dates and no definitive features are seen in cither case. 
We summarize the predicted fluxes and limiting sensitiv- 
ities for each low-redshift interloper case in Table 5. 

In the case of Abell 68, if the J-band line was Ha, 
for our adopted template H II region spectra, [Oil] and 
H/3 would have been recovered with > 45a in the LRIS 
spectrum. 

If the J-band feature was [Oil], then Lyman-a would 
have been seen with ease in the blue LRIS spectrum. 
While the resonant nature of the transition makes it dif- 
ficult to predict the Lyman-a flux , if even 0.1% of the 
photons can escape, a ~5cr detection would still have 
resulted. Since Lyman-a emis sion is often not see n in 
distant Lyman break galaxies (IShaplev et al.l 120031 ) . we 
can secure additional constraints from the non-detection 
of H/3 and the two [OIII] lines in the H-band spectrum. 
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Although some of the relevant regions are affected by 
OH emission, [OIII] 5007 A would be strong enough to 
be detected at >5<r under some cases. 

If the line were [OIII]A5007A , then H/3, [OIII]A4959A , 
and Ha would all fall in the spectral region covered by 
the H-band spectroscopy. H/3 and Ha would each be seen 
at >5 a for most of the range of emission line ratios. 

One seemingly unlikely hypothesis deserves attention. 
The candidate is quite close (1.8 arcsec) to a z=1.58 
galaxy with H/3 and [OIII] revealed in the J-band discov- 
ery spectrum (Richard et al. 2006, in preparation). If the 
candidate line were H/3, the rest-frame velocity separa- 
tion between the two sources would only be 400 km s -1 . 
The likelihood of an association thus seems high. How- 
ever, it turns out the combined J+H spectra of the two 
sources are very different. Ha is seen from the foreground 
galaxy in the H-band but the same line can be rejected for 
the candidate with 15a confidence for a Case B Balmer 
line ratio. Likewise the limiting sensitivity at the loca- 
tion of [OIII] emission line for the foreground hypothesis 
is sufficiently deep such that we would have expected to 
see a line for nearly all reasonable line ratio assumptions. 
In combination, therefore, the non-detection of Ha and 
[OIII] A 5007 A strongly rule out any association of the 
candidate with the z=1.58 galaxy. 

Similarly, in the case of Abell 2219 cl if the J-band 
line were Ha, [Oil], H/3 and [O III] 4959A would have 
been seen in the LRIS spectra at 23-85, 8.5, and 6.4-43 
a respectively. 

If the line was H/3, Ha would have been detected at 
22(7 in the H-band spectrum, and [OIII] A4959A and 
A5007A at 3.5-24 and 21-190 a respectively in the J- 
band. 

If the line were [Oil], then, as before strong Lyman- 
a would be expected in the blue LRIS spectrum, even 
if only a small fraction of photons escaped. And 
[OIII]A4959A and A5007A would have been seen at up 
to 4.5-8.0 and 4.8-ller respectively in the H-band spec- 
trum. 

Finally, if the line were [OIII] A5007A , then H/3 would 
be detected at 1.7— 15cx, depending on the oxygen abun- 
dance, and Ha would been seen in a clean area of the 
H-band with 4.8-43<r. 

All cases described above assume an HII region tem- 
plate spectrum. Assuming the emission lines are in- 
stead generated by a harder ionization source such as an 
AGN, the CIV AA 1548, 1550 line may be prominent in 
the rest-frame UV spectrum. If Lyman-a is suppressed 
via resonant scattering and absorption by dust, narrow- 
line AGN can exhibit line ratios of 1< CIV/Lya < 1.4 
(iDev et all 119951 : iDe Breuck et al.l 120011 : iDawson et all 
I2003J). If these sources arc indeed Lyman-a emitters at 
z ~ 10, then the CIV line may also be observable in 
the H-band. For both Abell 68 cl and Abell 2219 cl, 
CIV would lie in clean regions (1.60 /im and 1.55 fim, 
respectively); however, neither emission line is seen, sug- 
gesting that if the candidates are at z ~ 10, they have 
CIV/Lya < 1. 

5.6. Other Candidates 

Comprehensive optical and H-band spectroscopy is 
also available for Abell 2219 c2 and in fact the constraints 
we discuss above for Abell 2219 cl are just as convincing 



here making it a very promising source. However, we 
prefer not to elevate it into the status of the two most 
promising sources because of the potential broad-band 
ACS detection (§5.2). In this sense it represents a source 
of status intermediate between our top two candidates 
and those discussed below for which, due to weather and 
other problems, our spectroscopy is currently incomplete. 

For Al689cl, c2 and c3, we have no optical or H-band 
spectroscopy and so major interloper hypotheses such as 
Ha and [O II] cannot yet be excluded. 

In summary, we can conclude that at least two (and 
quite possibly three) of our six candidates cannot be ac- 
counted for as lines from lower redshift emitters. More- 
over, there is no reason to denigrate the status of the 
other 3 candidates either since, so far, every additional 
spectrum we have taken has continued to support the 
Lyman-a hypothesis. Our low contamination ratio is 
not that surprising when one considers all of the con- 
sidered interloper hypotheses would still represent lensed 
systems and thus, if highly magnified, not be seen at the 
location of the z ~10 critical line. 

On the other hand, it is clear that our tests make nec- 
essary assumptions about the line ratios expected from 
0.85 < z <2.26 galaxies. It is conceivable that these as- 
sumptions are invalid and that the sources have anoma- 
lous line ratios (although it would seem surprising that 
this is the case in several systems). Given the wavelength 
ranges accessible to ground-based spectrographs, we be- 
lieve we have undertaken all possible tests. 

6. IMPLICATIONS 

Noting that at least two of our six candidates may in- 
deed lie at z ~10, we now address the goal of the survey: 
what is the volume density of high z emitters? How- 
ever, noting their identification could still be considered 
tentative, we will also explore the implications if none 
of our candidates is real. After making simple assump- 
tions about the effect of dust and resonant scattering on 
escaping Lyman-a photons, we will use the derived abun- 
dance to constrain the star formation rate density. We 
then attempt to understand whether the derived abun- 
dance represents a significant contribution to the ionizing 
photon budget at z ~ 10. 

6.1. Volume Density of z>8 Lyman-a Emitters 

In a blank-field survey, the surface density of sources 
brighter than a given limiting line flux can be computed 
by dividing the number of detected sources by the area 
observed. The calculation is more complex in this case 
because of the varying amplification of the source flux 
over the survey volume and the distortion, by lensing, of 
the survey area. 

The differential comoving volume element at (tt,z) 
along one of the NIRSPEC slit positions is 



dV c (ft,z) 



1 



M(n,z) 



dl c (z) 



dz 



[D 2 c (z)dSl]. (3) 



where dl c is the comoving length of the volume element 
along the line of sight. We also compute the uncertainty 
associated with each volume element due to the error in 
the magnification (§3.4). 

The total survey volume over which sources with 
Lyman-a line luminosity greater than L is the integral 
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over all volume elements with limiting Lyman-a lumi- 
nosity below L, 

V C (L)= [ [ dV c (n, z) H[L - L nm (n,z)], (4) 

where H[ L-L\ lm (n, z) ] is a step function equal to 1 when 
L> Lii m (l~2, z) and equal to when L< Ln m (Cl, z). 

The uncertainty in the magnification introduces error 
into the derived comoving volume. We quantify this error 
by computing the distribution of comoving volume for all 
acceptable cluster mass models (see § 3.4). The range of 
comoving volume that spans ±34% of the total samples 
is taken to be the uncertainty in the volume. 

The survey volume and associated errors are plotted 
as a function of Lyman-a luminosity in Figure 10. The 
comoving volume decreases from ~30 Mpc 3 at L=3x 10 42 
erg s" 1 to ~3 Mpc 3 at L=3xl0 41 erg s -1 . The 'one- 
sigma' error bars in the volume are small enough that 
they do not extend outside the datapoints. In Figure 11, 
we plot the number density of z ~9 Lyman-a emitters 
as a function of Lyman-a luminosity assuming that (i) 
none of our candidates is at z ~8-10 (left) and (ii) two of 
our candidates are Lyman-a emitters (middle) and (iii) 
all six of our candidates are real (right). 

In the case where two of the candidates are real, the 
Poisson errors are 100%. This arises because galaxies 
need not only be brighter than the survey limiting lumi- 
nosity but must also lie within the area covered by the 
survey. In lensing experiments, in order for inclusion as- 
sociated with some limiting luminosity, the 5a limiting 
luminosity at the location of the object must be lower 
than the limiting luminosity of the survey. The net re- 
sult is that objects with fluxes close to the sensitivity 
limits only are included in "surveys" with a small range 
of limiting luminosities. If two candidates are real, as 
considered here, the luminosity ranges over which the 
objects would be detected do not overlap; hence no bin 
contains more than one object. 

To place the inferred abundances in context, we com- 
pare our results with other surveys that have been con- 
ducted, particularly at z~6 where the luminosity func- 
tion of Lyman-a emitters is more firmly established. 

The first observational constraint on the abun- 
dance of Lyman-a emit ters at z>8 was presented by 
( Willi s fe Courbir] l2005f ) . These authors conducted a 
narrow-band imaging survey in the J-band toward the 
Hubble Deep Field South and found no z = 8.8 sources 
with Lyman-a luminosity greater than 10 42 ' 5 erg s" 1 over 
a volume o f appro ximately 990 Mpc -3 . More recently, 
ICuby et "all (|2006| ) reported a null detection of z = 8.8 
Lyman-a emitters brighter than 10 43 erg s _1 over a co- 
moving volume of ~ 4200 MjPc 3 . Given the large error 
bars, the upper limits from I Willis fc Courbinl {2005) and 
ICuby et all 42006) are consistent with the case in which 
only two of the candidates are real. If all six of the can- 
didates are Lyman-a emitters, it becomes slightly more 
difficult to explain the lack of Lyman-a emitters in the 
two narrowband surveys without resorting to luminosity 
dependent evolution. Hence, in this case, the observa- 
tions could suggest that the decline in luminous Lyman-a 
emitters relative to fainter sources observed at z ~ 6 — 7 
(IShimasaku et alj|2006t iKashikawa et al.ll2006t live et all 
|200jJ continues to z ~ 10. 



A second question is whether our density is higher 
or lower than the Lyman-a luminosity function which 
is now fairly w ell-established at z = 5.7, at least at 
the b right end (ISantos et al.ll2004t iMalhotra fe Rhoadsl 
[2003 IShimasaku et al.ll200l . It is worth remembering 
that despite the fact that over 100 Lyman-a emitters are 
now catalogued at this redshift, the degree of spectro- 
scopic confirmation is still quite limited. The extant data 
are well-fit by a Schechter function with (j>*—1.6x 10 -4 
Mpc -3 and L*=1.6x 10 43 erg s _1 where, for the purposes 
of discussion, we have adopted the steepest likely faint- 
end slope of a = —2 so as to provide the most optimisitc 
no-evolution scenario. Assuming, naively no evolution 
between z=5.7 and z=10, it is unlikely that we would 
find a source brighter than 10 40 erg s _1 . At face value, 
therefore, our detection of at least two promising cases 
would suggest an upwards evolution in the abundance of 
faint emitters with redshift. 

The likely evolution of the Lyman-a luminosity func- 
tion between z=5.7 and z~10 is a complicated func- 
tion of the star formation efficiency, stellar initial 
mass function, relative geometry and kinematics of HII 
and HI gas, and dust content of the host galaxy, as 
well as the ionization state of the IGM. If, for ex- 
ample, the dust content of galaxies decreases toward 
higher redshifts, the percentage of Lyman-a photons 
that escape the host galaxy may increase. On the 
other hand, as the IGM becomes increasingly neu- 
tral, observed Lyman-a count s are expected to de - 
cline (jMalhotra fe Rhoadsl l200l lHaiman fc Ced 120051 ) . 
although the presence of large HII regions around clus- 
tered star- forming sou rces may slightly offset the mag- 
nitude of this decline (jFurlanetto et al.1 12004). Clearly, 
it is very difficult to predict exactly how the Lyman-a 
luminosity function will evolve over this ~500 Myr time 
interval. 

Bearing these uncertainties in mind, our detection of 
two promising candidates clearly implies a steep faint- 
end slope and hence a large abundance of intrinisically 
faint sources. If more of our candidates are at z ~ 10, it 
would almost certainly imply that low-luminosity star- 
forming systems are more abundant at z — 9 than at 
z = 5.7. This may be a result of evolution in galaxy 
properties (as discussed), or the effect of a predominantly 
neutral IGM which increases star formation efficiency in 
low-mass halos due to the decrease in the Jeans mass 
(Bark ana fc Lo cb 2001). Determining which of these fac- 
tors is primarily responsible for the possible upturn in 
counts is difficult with the available data. 

If none of our candidates is real, a number of explana- 
tions are possible. First the faint-end slope may not be 
as steep as required for us to detect sources in our ob- 
served volume. There are a number of feedback processes 
(galactic winds, supernova explosiona, etc.) that could 
significantly reduce star formation efficiency in low-mass 
halos, thereby resulting in a relatively shallow faint-end 
slope. Alternatively, if the IGM is significantly neutral, 
it is possible that faint-end slope is very steep, but the 
neutral IGM is blocks our view of these systems. We de- 
fer a detailed discussion of the implications of our survey 
in terms of what is now known at z ~5-6 to a later paper 
(Stark, Loeb & Ellis, 2007). 



A Keck Survey for Gravitationally-Lensed Lyman-a Emitters in at 8.5<z<10 



11 



6.2. Contribution of Low Luminosity Galaxies to 
Reionization 

Several methods have been adopted for considering the 
contribution of star-forming sources to cosmic reioniza- 
tion. One approach is to compare the inferred ionizing 
photon production rate from star-forming galaxies of all 
luminosities to that required to reionize neutral hydro - 
gen in the IGM (jMadau et al .11 1 9 9 9t IStTavelli et alJl2004ft . 
Since the abundance of Lyman-a emitters at z ~ 10 is 
very uncertain, comparing the total ionizing photon out- 
put is highly dependent upon assumptions about the lu- 
minosity function. In this initial study, we ask a more 
basic question: do low luminosity Lyman-a emitters pro- 
vide a significant contribution to reionization? Following 
this approach, we therefore only tabulate the ionizing 
photon production arising from low luminosity (i.e. 10 41 
- 10 42 erg s _1 ) sources that our survey is designed to 
probe. 

The contribution of low luminosity star-forming 
sources to reionization can be parameterized in simple 
terms by comparing their ionizing photon production 
rate to that required for reionization: 

-o ^B U f cscr y m* \-i 

ngal ~ Z \10J\0m) Vo.lM yr-iJ 

assuming the space density of sources is constant during 
the reionization epoch. Here, B is the number of ion- 
izing photons required to keep a single hydrogen atom 
ionized, njj is the comoving number density of hydrogen 
at the redshift of interest, At is the period over which 
reionization occurs, f esc is the escape fraction of ionizing 
photons, M* is the time-averaged star formation rate of 
galaxies in units of M Q yr _1 . 

In both cases, we fix the comoving number density of 
hydrogen at nn=1.9 x 1CP 7 cm -3 , appropriate for the 
comoving hydrogen number density adopting the best-fi t 
WMAP cosmological parameters (|Spergel et alj [2006) . 
We have also assumed an IMF with Salpeter slope with 
stellar masses ranging from 1 M Q to 100 Mm and metal- 
licityofof Z=0.001 (l/20th solar) (|Schaererll2003f) . If the 
IMF is weighted more toward massive stars or is substan- 
tially more metal-poor, the ionizing photon output for a 
given star formation rate will be greater, and less galaxies 
will be required for reionization. 

The escape fraction of ionizing photons is not well- 
constrained observationally; we allow this quantity 
to vary between 02 and 0.5 dSteidel et all 120011 : 
iFernandez-Soto et alJl2003t IShaplev et alJl2006f T We as- 
sume the escape fraction is constant with luminosity; 
however, models of high-redshift galaxies suggest that 
it should increase to ward the low-luminos ity end of the 
luminosity function (jRicotti fc S hrill 2000) . If this is the 
case, then contribution of low-mass galaxies to reioniza- 
tion will be even greater than what is shown below. 

The recombination rate increases with the local hy- 
drogen number density, hence the number of ionizations 
required per atom for reionization depends strongly on 

<n 2 > 

the dumpiness of the IGM (C= <w ' ? > a ). Simulations 
gener ally predict dumpiness factors of C~30 (|Gnedinl 
2000) for z^6-7. Following the interpolation formulae 



in (|Stiavelli et a l. 2004), this results in ~5-10 ionizations 
per hydrogen atom depending on the timescale of reion- 
ization. 

Substituting the range of parameter values into equa- 
tion 5, we find that for galaxies with star formation rates 
of 0.1 Mq yr _1 to dominate the reionization process, they 
should have comoving number densities in excess of 0.1 
- 4 Mpc~ 3 (Figure 12). If two or more of the candidates 
are at z~10, then it would appear that low luminosity 
galaxies contribute significantly to the reionization pro- 
cess. 

While there is currently not enough data at z ~ 
10 to estimate the integrated ionizing photon output 
from Lyman-a emitters, we can attempt to constrain 
this quantity by assuming that the luminosity function 
evolves only in normalization (hence maintaining the 
same shape) between z = 6.5 and z — 10. In the most 
optimistic case where six candidates are at high-redshift, 
we take the characteristic lumi nosity, L* derived f or the 
z = 6.5 Lyman-a emitters by [Kashika wa et all (2006) 
with a diverging faint-end slope of a = —2. Integrating 
the luminosity function down to the survey limit (~0.01 
Mq yr _1 ), we find that the total ionizing photon produc- 
tion rate would be sufficient for reionization for nearly 
the full range of assumptions considered in Figure 12. 

7. SUMMARY 

Low luminosity star-forming galaxies are often pre- 
dicted to produce a large supply of ionizing phot ons dur- 
ing the reionization era (|Barkana fc Loeb1l2001f ). How- 
ever, feedback processes from supernova explosions or 
galaxy winds may significantly decrease star formation 
efficiency in the low-mass dark matter halos that are 
thought to house the majority of low-luminosity star 
forming sources. Tentative observational evidence hints 
at a very large abundance of low luminosit y sources 
(|Yan fc Windhorst! 120041 : iBouwens et~aT1 I2006D suggest- 
ing that these feedback processes may not significantly 
affect star formation in low-mass halos. Until the con- 
struction of JWST and 20-40 meter class ground-based 
telescopes, gravitational lensing surveys offer one of the 
most efficient means of spectroscopically verifying the 
presence of an abundant population of low-luminosity 
(< 10 43 ergs _1 ) sources. With this as our goal, we 
have conducted a spectroscopic survey for gravitation- 
ally lensed Lyman-a emitters over the redshift interval 
8.5< z <10.4. We summarize our results below. 

1. The critical lines of 9 well- understood lensing 
clusters were scanned spectroscopically in the J-band 
for lensed Lyman-a emission from galaxies at z=8.5 to 
10.4. The new observations are sensitive to star-forming 
sources with intrinsic (unlensed) Lyman-a luminosities 
of 10 415 erg s~ x over a significant fraction of the survey 
area. This limit is over an order of magnitude fainter 
than the corresponding sensitivity limits of conventional 
surveys for Lyman-a emitters at z > 6. 

2. Six promising candidate Lyman-a emitters were 
identified whose reality we justify noting the possibility 
of instrumental and detector artefacts. The emission fea- 
tures have fluxes that range from 2-5xl0 -17 erg cm~ 2 
s _1 and significances between 5-8<r. The Lyman-a lu- 
minosities of the sources range from 10 41 to 5xl0 42 erg 

3. We attempt to verify that our lines arise from z ~10 
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Lyman avia several tests. No convincing broadband op- 
tical counterparts were seen in deep HST images (R~27) 
except in one marginal case. For 3 out of the 6 cases, 
we conduct exhaustive additional spectroscopy to explore 
the likelihood that low-redshift interloper lines may ex- 
plain the J-band emission. In no case do we find such an 
example, and, overall we conclude that at least 2 of our 
6 candidates probably lie at z ~10. We briefly discuss 
the considerable challenges of making further progress in 
confirming any or all of our candidates given the limited 
rest-wavelength range available with large ground-based 
telescopes. 

4. Assuming two or more of our Lyman-a emitter can- 
didates are real, then the cumulative abundance of low 
luminosity galaxies (defined as those with L>10 41 ' 5 erg 
s _1 ) is at least 0.3 Mpc -3 . Such a large abundance of low 
luminosity Lyman-a emitters supports the contention of 
a steep faint end slope for the star-forming luminosity 
function at z~10. 

5. Assuming a single population of galaxies dominates 
the reionization process the number density of galaxies 
forming stars at a rate of 0.1 M Q yr _1 must exceed 0.1-4 
Mpc~ 3 to reionize neutral hydrogen. Our inferred abun- 
dance is consistent with being within this range indicat- 
ing that intrinsically faint Lyman-a emitters may well 
produce a substantial supply of ionizing photons during 
the reionization era. 

6. However, if none of our candidates are at high- 
redshift, then there are several likely explanations. If the 
IGM is significantly neutral, low luminosity star form- 
ing systems may be very abundant but may not be ob- 
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APPENDIX 

The observing sequence for narrowband surveys often involves alternating between two narrowband filters (e.g. A 
and B) with slightly different central wavelengths. In each filter, the observations are divided into subexposures with 
integration times that are comparable to the variability timescale of the sky background (typically 5-10 minutes). 
A composite image is made of all of the subexposures in each filter, resulting in two "subsurveys" slightly offset in 
redshift space. Subtracting the two composite images removes continuum soures, thereby allowing line emitters to be 
identified. Alternatively , instead of a second narrowband filter, a broadband filter can be used to identify line emitters 
( Willi s~fc Courbinl feOOS). Using either approach, if integration time tint is devoted to a narrowband survey, then the 
total integration devoted to each individual "subsurvey" is only (t; n t/2) in length. In the following calculations, we 
assume a total integration time of 60 hours. 

The field of view of the survey camera is one of the key parameters in determining the survey sensitivity. For a single 
exposure, the field of view is limited by the size of the detector and the plate scale of the detector. The plate scale is 
typically set such that each pixel samples 1/4 of a spatial resolution element. Assuming seeing-limited observations 
with a 0'.'6 seeing disk, this corresponds to 0.1-0.2 arcseconds per pixel. Until very recently, only 1024x1024 pixel 
detectors were available on 8-10 meter telescopes in the near-infrared; given typical plate scales, this allows fields of 
view between 3 and 12 arcmin 2 . Recently, larger format cameras have become available (e.g. MOIRCS on Subaru) 
allowing fields of view of ~30 arcmin 2 to be observed in single exposures. In the signal to noise calculations that 
follow, we consider both a detector with a field of view of 4 arcmin 2 (similar to that of NIRI on Gemini) as well as 
one with a field of view of 30 arcmin 2 . 

In the near-infrared, the noise is dominated by bright atmospheric emission lines from OH molecules; hence the 
background in a given exposure is dependent on the central wavelength and bandwidth of the narrowband filter. The 
central wavelength is chosen to lie in >50A gaps between bright OH lines, where the atmospheric transmission is 
near 100%. We assume a central wavelength of 11905 A for our predictions. The signal to noise in a single exposure 
is minimized when the filter bandwidth is matched to the FWHM of the emission line, which is typically 5-10 A for 
Lyman-a emitters at high redshift (Shimasaku et al. 2006); hence, we adopt 10A as the FWHM for our calculations. 
We also consider a filter with a FWHM of 100 A; while the noise in a single exposure is larger for this filter, the 
redshift range covered is ten times larger (Az=0.1 vs. Az=0.01). 

The received signal from the science source is estimated as follows, 

S source = J f X „r(X) dX A tMin ^ QEVa UU nt /2) (A-l) 

where i\, source is the spectrum of the source in units of photons s^ 1 cm -2 A -1 , t(A) is the filter transmission profile 
(assumed gaussian with FWHM defined above), rj te i is the telescope optics throughput, r\i ns t is the throughput of the 
instrument, t\qe is the quantum efficiency of the detector, rj atm is the atmospheric transmision, A te i is the collecting 
area of the telescope. As discussed in §5.1, we expect the line profile of high-redshift Lyman-a emitters to be asymmetric 
with the blue side of the line absorbed. Along these lines, we assume the line profile is a half-Gaussian with no flux 
shortward of the central wavelength. While this profile is perhaps oversimplified, the S/N calculations do not change 
signficantly if the profile is altered. We adopt throughput and quantum efficiency values that are consistent with 
those expected for DAzLE (?7 fe ;=0.6, 7 7i ras t=0-5, t]qe=0.7), a near infrared narrowband imager designed to find z > 7 
L yman-q emit ters iHorton et al.l (|2004h . The atmospheric transmission, T] atm , is computed following models presented 
m iLordl (I1992T) assuming an air mass of 1.5 and a 1.0 mm water vapour column 5 and is near unity for the OH window 
considered. We assume a telescope collecting area of 7.9 x 10 5 cm 2 corresponding to the size of Keck. The photometric 
aperture is matched to the area subtended under reasonably good seeing conditions (~ 076 diameter seeing disk). The 
integration time is only (ti nt /2) because the total integration time is split between two different narrowband filters, as 
described above. 

The noise in an exposure of time tint is a function of the integrated source counts, sky background, detector dark 
current, and detector read noise. Formally, the noise is defined as 

N = yjN sk y + N source + N DC + RN 2 (A-2) 

where N s k y , N source , N^c a re the number of electrons from the sky background, science source, and dark current, and 
RN is the read noise. 

We compute the expected background count rate utilizing the Gemini model for the near-infrared sky background 
spectrum 6 . The model inludes the contribution from atmospheric OH emission lines, zodiacal emission (from a 5800 
K blackbody), and thermal emission from the atmosphere (from a 250 K blackbody). The total counts from the sky 
background in a photometric aperture of area & p hot in the composite image, is given by 

N sky = J h , sky r(X) dX rj tel rj inst rj QE A tel a phot (t int /2) (A-3) 

where f\. s ky is the sky background in photons s _1 cm~ 2 A" 1 arcsec -2 , r is the filter transmission profile (assumed 
gaussian with FWHM defined above), & p hot has units of arcsec 2 , and all other parameters are the same as defined 
above. 



5 available from http : //www. gemini . edu/sciops/0bsProcess/obsConstraints/atm-models/trans_10_15 . dat 

6 available from http : //www. gemini . edu/sciops/0bsProcess/obsConstraints/atm-models/nearIR_skybg_16_15 .dat 
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TABLE 1 

Clusters surveyed 



Cluster 



Rcdshift 



RA a 



Dec b 



Lens Model Reference 



CI 0024 
Abell 68 
Abell 370 
MS 0451 
Abell 963 
Abell 1689 
Abell 2218 
Abell 2219 
Abell 2390 



0.39 
0.255 
0.375 
0.55 
0.206 
0.183 
0.176 
0.226 
0.228 



00 26 35.5 
00 37 06.8 
02 39 53.1 
04 54 10.6 
10 17 03.7 
13 11 29.4 
16 35 49.3 
16 40 19.8 
21 53 36.9 



+17 09 50.7 
+09 09 23.4 
-01 34 54.8 
-03 00 50.7 
+39 02 49.2 
-01 20 28.7 
+66 12 43.5 
+46 42 41.9 
+17 41 43.4 




Referen ces. — fltlKneib et ail ll20u3h: (2) ISmith et ail (2005ft (3) IBezecourt et ail (l99Sft : (4) IBorvs et all (200lft . (5) Limousin et al. 
in prep; (6l lKneib et a .1 (|20041 1: (7) ISwinbank et ail (2006ft 
a units of HH MM SS b units of +DD MM 

The total counts associated with the dark current is given by 



where n p i X is the number of pixels in the photometric aperture area and dc is the dark current in units of electrons 
per second per pixel. The number of pixels in the aperture is derived assuming a plate scale of 0'.'15. The typical dark 
current for near-IR imagers on 8-10 meter telescopes is 0.25 electrons s" 1 pixel -1 ; hence we adopt this value in our 
calculations. 

The read noise contribution to the total noise is computed assuming Fowler sampling is used, which reduces the 
read noise by a factor of y/4/n rea d where n rea( j is the number of reads at the beginning and end of each exposure. The 
total read noise in an aperture of size n p j x in the composite images is thus given by 



where rn is the detector read noise (for single readout) in units of electrons per p ixel and n exv , is th e number of 
exposures in the composite frame. We adopt a read noise of 10 electrons per pixel (|Horton et alj|2004l ) for a single 
readout and assume the detector is read out 16 times per exposure. 

The 5<7 limiting flux is defined as the line flux for which the signal from the science source is a factor of five greater 
than the noise. For the instrument properties assumed, the 5a limiting flux in a 10 hour narrowband survey is 3- 
9xl0 -18 erg cm~ 2 s _1 for a 10A and 100A narrowband filter, respectively. At z=9, this corresponds to a Lyman-a 
luminosity of 0.6-2xl0 42 erg s _1 assuming a magnification of x5. The survey volume is computed following equations 
3 and 4 presented in §5.1.2 and assuming a magnification distribution for one of the most well-studied clusters, Abell 



We also compute the efficiency of a conventional narrowband survey. The survey parameters are identical to those 
assumed for the lensing survey. The only difference is that in computing the limiting luminosity and survey volume, 
we do not fold in the magnification due to the lensing cluster. 

The parameters chosen for the spectroscopic survey correspond to the properties of NIRSPEC on Keck. The limiting 
sensitivity in a spectroscopic survey is computed in the same way as for the narrowband survey. The spectroscopic 
survey strategy is described in §3.2. We divide the total integration time into many 1.5 hour observations; each 1.5 
hour integration is further divided into nine 10 minute exposures. One key difference between the spectroscopic and 
narrowband survey is that the entire integration time is devoted to a single survey; hence, the (tj nt /2) in the equations 
above become U n t for spectroscopic surveys. 

Most of the parameters used in the spectroscopic S/N computation are identical to those used above. The read 
noise in each ten minute exposure is 25 electrons pixel -1 assuming 16 total reads; while slightly higher than in the 
narrowband survey described above, the noise is still sky-background dominated. 

The median 5a limiting flux for the spectroscopic survey is 8xl0 -18 erg cm" 2 s _1 . We assume the longslit is 
placed along the critical lines of the lensing clusters listed in Table 1 , which results in a significantly higher median 
magnification than in the narrowband survey. For a magnification of x25, this corresponds to a limiting Lyman-a 
luminosity of 4xl0 41 erg s _1 . 



Ndc = n ptx dc(U nt /2) 



(A-4) 




(A-5) 



1689. 
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TABLE 2 
NIRSPEC Survey Observations 



Date 


Cluster 


Position Anglc a 


Integration timc b 


Aug 2004 


Abcll 2390 1 


117 


12.0 




Abell 2219 1 


119.2 


9.6 




Abell 68 1 


330. 


13.2 


Jan 2005 


MS 0451 1 


120.9 


11.4 




Abell 963 1 


177.2 


10.8 




Abell 963 2 


175.0 


5.4 




Abell 1689 1 


43 


5.4 


June 2005 


Abell 1689 1 


43 


3.6 




Abell 1689 2 


205 


10.8 




Abell 1689 3 


139.1 


10.8 




Abell 2218 1 


153 


10.8 




Abell 2218 2 


134 


5.4 




Abell 2219 1 


119.2 


5.4 




Abell 2219 2 


160 


16.2 




Abell 2390 2 


109 


10.2 




Abell 2390 3 


20.0 


6 


Oct 2005 


Abell 68 2 


300 


10.8 




CI 0024 1 


139.7 


5.4 




CI 0024 2 


160.2 


4.8 




MS 0451 2 


100.1 


7.2 




Abell 2390 3 


20.0 


5.4 




Abell 2390 4 


90.0 


4.8 




Abcll 2219 2 


160 


2.4 



a in degrees North through East b in ksec 



TABLE 3 
NIRSPEC Lyman a Candidates 



Candidate 


RA(J2000) 


DEC(J2000) 


flux 

(erg s" 1 cm" 


- 2 ) 


A 
(/urn) 


z Lya 


LoglO L Lya, 

(erg s" 1 ) 


Logio Mag 


A68cl 


00 37 06.10 


+09 09 18.7 


2.2+0.3x10" 


17 


1.254 


9.32 


41 2+ ai 
4i - z -0.7 


n 9 +0.7 
z ' z -0.1 


A1689cl 


13 11 31.85 


-01 20 40.0 


2.3+0.3x10" 


17 


1.366 


10.23 


42 n+0.02 


i 9+O.OI 


A1689c2 


13 11 31.26 


-01 20 29.8 


4.0+0.7x10" 


17 


1.173 


8.65 


42 7 +0.07 
^■'-0.08 


n a+001 
U - 8 -0.03 


A1689c3 


13 11 28.43 


-01 19 44.5 


4.1+0.6x10" 


17 


1.259 


9.35 


42 -,+0.06 
4z -0.08 


i c+0.04 
L -°-0.03 


A2219cl 


16 40 23.64 


+46 42 26.5 


4.8+0.6x10" 


17 


1.215 


8.99 


41 o+0.09 
41.8_ 07 


1 - a -0.07 


A2219c2 


16 40 23.04 


+46 42 43.2 


2.1+0.3x10" 


17 


1.209 


8.94 


42 n +0.06 
4z - u -0.09 


-, 0+O.O6 
L -°-0.03 



TABLE 4 

Wavelengths of Additional Emission Lines for Low-z Scenarios 



Line 


z 


A Lya 


A [on] 


A H/3 


A [oiii] 


A Ha 






w 


(/urn) 


(H 


(/im) 


H 


Ha 


0.85 


0.2253 


0.6905 a 


0.9007 


0.9188/0.9277 


1.2160 


[OKI] 


1.43 b 


0.2953 


0.9051 


1.1805 


1.2043/1.2160 


1.5939 




1.45 c 


0.2982 


0.9139 


1.1920 


1.2160/1.2278 


1.6093 


H/3 


1.50 


0.3042 


0.9323 


1.2160 


1.2405/1.2525 


1.6418 


[OH] 


2.26 


0.3967 


1.2160 


1.5860 


1.6180/1.6336 


2.1413 


Lya 


9.0 


1.2160 


3.7270 


4.8610 


4.9590/5.0070 


6.5630 



a Throughout the table, we take the luminosity- weighted value of [Oil], A=3727 A, for conciseness. b A=5007 A c A=4959 A 
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TABLE 5 

Constraints on Low-redsrtft Interpretations of Candidates 



Line 


z 




f [OH] 


f H/3 


f [om]i/ f [oiii] 2 


f Ha 






(10~ 17 erg cm~ 2 s 


i x ) (10 17 erg cm" 


- 2 s- 1 ) (10- 17 erg cm" 2 s" 1 ) 


(lO" 17 erg cm" 2 s" 1 ) 


(10~ 17 erg cm- 2 s- 1 ) 


Abcll 68 cl 


Ha 


0.91 


N/A 


0.63-2.4 (0.07) 


0.79 (0.07) 


N/A / N/A 


N/A 


[OKI] 2 


1.51 


N/A 


N/A 


0.31-2.8 (1) 


0.73 (0.9)/ NA 


0.4-13(2) 


[OIII]i 


1.53 


N/A 


N/A 


1.1-7.7 (6) 


NA / 6.9 (5) 


3.1-21 (2) 


H/3 


1.58 


N/A 


N/A 


NA 


0.66-4.4 (5) / 1.8-15 (3) 


6.2 (2) 


[Oil] 


2.37 


31 (0.02) 


N/A 


1.4 (5) 


0.22-7.3 (7) / 2.2-5.1 (5) 


N/A 


Abell 1689 cl 


Ha 


1.06 


N/A 


N/A 


N/A 


N/A / N/A 


N/A 


[OKI] 2 


1.71 


N/A 


N/A 


0.33-2.9(6) 


0.77(1) / N/A 


N/A 


[OIII]i 


1.72 


N/A 


N/A 


1.2-7.7(0.9) 


N/A / 6.9(4) 


N/A 


H/3 


1.79 


N/A 


N/A 


N/A 


0.69-4.6(20) / 1.8-16(100) 


N/A 


[OH] 


2.64 


N/A 


N/A 


N/A 


N/A / N/A 


N/A 


Abell 1689 c2 


Ha 


0.78 


N/A 


N/A 


N/A 


N/A / N/A 


N/A 


[OKI] 2 


1.33 


N/A 


N/A 


N/A 


1.3 (5) / N/A 


N/A 


[OIII]i 


1.36 


N/A 


N/A 


N/A 


N/A / 12(3) 


N/A 


H/3 


1.40 


N/A 


N/A 


N/A 


1.2-8.0(1) / 3.2-38(4) 


N/A 


[OH] 


2.13 


N/A 


N/A 


N/A 


N/A / N/A 


N/A 


Abell 1689 c3 


Ha 


0.94 


N/A 


N/A 


N/A 


N/A / N/A 


N/A 


[OIII] 2 


1.54 


N/A 


N/A 


0.59-5.1 (8) 


1.4 (2) / N/A 


N/A 


[OIII]i 


1.57 


N/A 


N/A 


6.2-11 (2) 


N/A / 12 (5) 


N/A 


H/3 


1.62 


N/A 


N/A 


N/A 


1.2-8.2 (3) / 3.3-29 (1) 


N/A 


[OH] 


2.42 


N/A 


N/A 


N/A 


N/A / N/A 


N/A 


Abell 2219 cl 


Ha 


0.85 


N/A 


1.4-5.1 (0.3) 


1.7 (1) 


0.51-3.4 (0.4) / N/A 


N/A 


[OIII] 2 


1.43 


N/A 


N/A 


0.69-6.0(2) 


1.6 (5) / N/A 


1.9-17(2) 


[OIII]i 


1.45 


N/A 


N/A 


2.4-16(0.9) 


N/A / 14(1) 


6.7-45 (7) 


H/3 


1.50 


N/A 


N/A 


N/A 


1.4-96 (2) / 3.8-34(0.9) 


13 (3) 


[OH] 


2.26 


35-130 (0.4) 


N/A 


1.6-6.0 (8) 


1.8-3.2 (2) / 4.8-11 (5) 


N/A 


Abell 2219 c2 


Ha 


0.84 


N/A 


0.6-2.3 (0.4) 


0.75 (1) 


0.23-1.5(0.4) / 0.60-5.3 (0.5) 


N/A 


[OIII] 2 


1.41 


N/A 


N/A 


0.30-2.7 (2) 


0.7 (3) / N/A 


0.84-7.3 (10) 


[OIII]i 


1.44 


N/A 


N/A 


0.40-13 (1) 


N/A / 6.3 (2) 


0.4-13 (1) 


H/3 


1.49 


N/A 


N/A 


N/A 


0.63-4.2(3) / 1.7-15(0.8) 


5.9(4) 


[OH] 


2.24 


15-58 (0.4) 


N/A 


0.70-2.6 (2) 


0.79-2.4(4) / 2.1-4.9 (5) 


N/A 


Note. 


— The fluxes arc predictions given the flux of the J-band feature and typical flux ratios for HII regions. The value 


in paranthesis is 


the limiting line flux at the corresponding wavelength (see t 


;3.2 for details on how the limiting flux is computed). [OIII]i is A4959Aand [OIII]2 is 



A5007A 
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Fig. 1. — Comparison of different survey strategies for locating Lyman-a emitters with current instrumentation in seeing-limitcd 
conditions. Solid circles represent the minimum source abundance that a 60 hour spectroscopic lensing observing campaign (over 10 clusters) 
could constrain (at 5cr); while open circles correspond to the minimum source abundance that a narrowband imaging lensing survey could 
constrain in the same exposure time. Both survey techniques probe orders of magnitude fainter than a conventional narrowband survey 
(open square), albeit over much smaller volumes. At luminosities to the left of the dashed vertical line, the spectroscopic approach covers 
a larger volume, and thus it is significantly more efficient than narrowband imaging in constraining the abundance of the faintest (< 1 Mq 

yr _1 ) Lyman-o emitters. 
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Fig. 2. — Survey clusters, with survey area and lensing critical curves. For each cluster, superposed on the HST /WFPC2 image are 
slit positions observed with NIRSPEC and the critical lines for a source at z = 8 (dotted lines). The solid curves bound areas where the 
magnification provided to z = 8 background sources exceeds a factor of 10. From left to right and top to bottom are Abell 68, Abell 370, 
Abell 963, Abell 1689, Abell 2218, Abell 2219, Abell 2390, CL 0024, MS 0451. 
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Fig. 3. — Limiting flux for detection of 5cr Lyman-a emission line for a typical NIRSPEC J band exposure. The top axis denotes 
the redshift corresponding to observed Lyman-a at the wavelength along the bottom axis. OH bands lie throughout the J-band 
spectrum, signficantly increasing the limiting flux at those wavelengths. The widest band occurs in the middle of the 
spectrum at ~ 1.26 — 1.28/im. 



A Keck Survey for Gravitationally-Lensed Lyman-a Emitters in at 8.5<z<10 21 

10000 i i 




10 20 30 40 50 

Distance along slit [arcsec] 



Fig. 4. — Lensing magnification provided to background sources as a function of NIRSPEC slit position and redshift. This figure shows 
the amplification at one of the slit positions along the critical line of Abell 2219. The solid line corresponds to the amplification provided 
to background sources at z=9 whereas the dotted and dashed line denotes that provided to sources at z=8.5 and 10.0, respectively. The 
error in lensing magnification is plotted as a function of NIRSPEC slit position for sources at z=9. The error does not vary significantly 
as a function of source redshift. Even allowing for uncertainty in the mass model, the magnification is still very high across the entire slit. 
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Fig. 5. — Distribution of lensing magnfication across our total survey area. Th e median ma gnification across our survey area is ~20, 
about a factor of 2 larger than in a similar optical spectroscopic lensing survey {Santos et al. 2004) (see §3.4 for explanation). A 
non-negligible fraction of each slit is magnified by larger than this factor. 
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Fig. 6. — Candidate Lyman-a emitters, each undetected in deep optical HST data. If Lyman-ct, the emission line rcdshifts range between 
z = 8.7 and z = 10.2. Line fluxes range between 2 — 5xl0 -17 erg cm -2 s — 1 implying intrinsic (unlensed) star formation rates of 0.2-5 Mq 
yr~ - 1 , for nominal source assumptions. 
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Fig. 7. — Normalized histogram of summed flux in randomly placed apertures in clean regions of sky-subtracted spectra. The dotted 
line represents value of a typical candidate Lyman-a emitter. The frequency of features with fluxes on par with the candidates is 0.02%, 
consistent with a ~ 4<r detection assuming a Gaussian distribution of noise. 



A Keck Survey for Gravitationally-Lensed Lyman-a Emitters in at 8.5<z<10 25 




Fig. 8. — Deep broadband images of the locations of the candidate J-band emission features. None of the candidates have definitive 
detections in the optical; although A2219c2 is very marginally detected in the ^' F ^c >0L p and ^F\60W bands. Three of the candidates are 
located within several arcseconds of brighter objects (A2219cl, A1689cl, and A1689c2); however in both cases, our registration is known 
to sufficient precision to rule out the possibility that the J-band emission features emanate from the bright objects. The lack of detection 
in very deep near-IR images suggests that if these candidates are at z ~ 9, then their rest-frame UV continuum emission is very faint. 
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Fig. 9. — Most likely emission features and corresponding redshifts of the J-band emission features. If the line lies in the rest-frame 
optical ([Oil], H/3, [OIII], Ho), then additional emission features should be present in optical or near-infrared spectroscopy. If, on the other 
hand, the line is Lyman-a at z=8.5-10, then additional powerful emission lines will have been redshiftcd into the mid-infrared, making 
confirmation very difficult. 
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Fig. 10. — Survey volume (comoving) sensitive to sources with Lyman-a luminosity L. The volumes observed are very small as a result 
of the small area of the spectroscopic slit and the decrease in survey area resulting from gravitational lensing. 
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Fig. 11. — Constraints on number density of 8.5<z<10.2 sources with Lyman-a luminosity brighter than L. Left: The 5<r upper limit to 
the abundance of low luminosity Lyman-a emitters is presented assuming none of our candidate Lyman-a emitters is real. Middle: The 
cumulative abundance of Lyman-a emitters if two of the candidates (A68cl and A2219cl) are real. Right: The cumulative abundance of 
Lyman-a emitters assuming each of the six candidates is a Lyman-a emitter. If even one candidate is real, the inferred abundance of low 
luminosity Lyman-a emitters at z ~ 9 would be very high. 
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Fig. 12. — Abundance of sources required for reionization assuming a single population of star-forming galaxies dominates the process 
(blue swath) compared the observed abundances (a) if all the candidates are real (open circles), (b) two of the candidates are real (solid 
circles), or (c) if none of the candidates are real (dotted line). The abundance constraints for case a and case b are identical for the two 
lowest luminosity bins plotted in the figure. The range for the abundance of sources required for reionization is determined assuming 
f c =0.5, At re i on =575 Myr, and B=10 for the lower locus of points and f c =0.02, At re i on =250 Myr, and B=5 for the upper locus of points. 
If two or more of the candidates are real, then it appears that intrinsically faint galaxies may play a dominant role in the reionization of 
neutral hydrogen. 



